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Abstract. Thermal fluctuations are known to play an important 
role in low-dimensional systems which may undergo incommensurate- 
commensurate or (for an accidentally commensurate wavevector) lock-in 
transitions. In particular, an intermediate floating phase with algebraically 
decaying correlations exists only in D = 2 dimensions, whereas in higher 
dimensions most features of the phase diagram are mean- field like. 

Here we will show, that the introduction of frozen-in disorder leads to 
strong fluctuation effects even in D < 4 dimensions. For commensurate 
wavevectors the lattice pinning potential dominates always over weak im- 
purity pinning if p < p c = 6/tt (D = 3), where p denotes the degeneracy 
of the commensurate phase. For larger p a disorder driven continuous 
transition between a long-range ordered locked-in phase and quasi-long- 
range ordered phase, dominated by impurity pinning, occurs. Critical ex- 
ponents of this transition, which is characterized by a zero temperature 
fixed point, are calculated within an expansion in 4 — D. The generaliza- 
tion to incommensurate wavevectors will be discussed. If the modulation 
in the quasi-long-range ordered phase has hexagonal symmetry, as e.g. for 
flux-line lattices, the algebraic decay is non-universal and depends on the 
Poisson ratio of the elastic constants. Weakly driven transport is domi- 
nated by thermally activated creep in both phases, but with different creep 
exponents. 

Incommensurate (I) phases appear in a large variety of systems (for a review see 
e.g. [1]). Examples are: 

(i) Charge density waves in quasi one and two dimensional conductors (e.g. in 
TTF-TCNQ, 2H-TaSe). 

(ii) Spin density waves (e.g. in CuGe0 3 ). 

(hi) Mass density waves in adsorbed monolayers (e.g. He, Kr on graphite or metal 
surfaces) or of reconstructed surfaces (e.g. of Mo). 

(iv) Polarization density waves in ferroelectrics with an incommensurate phase 
(e.g. in K 2 Se0 4 ). 

(v) Flux density waves in type II superconductors or Josephson junctions in an 
external field. 

Charge density waves are usually accompanied by a mass density wave like in 
superionic conductors or reconstructed metal surfaces. A common feature of the 
I phase is, that the wave vector q , which describes the spatial modulation of the 



density wave in the absence of any coupling to the lattices, varies continuously 
with the parameter of the system (e.g. temperature T, pressure p, chemical 
potential /i, magnetic field H etc.). If 27r/go is close to a multiple of the spacing 
of the underlying crystal lattice, i.e. if |g/p — qo| = \S\ < 5 C , commensurability 
effects may become important. Here g denotes a reciprocal lattice vector of the 
crystal and p is an integer. The modulation then may become commensurate (C) 
with the crystal lattice. The most striking effect of the C-phase is the existence of 
a gap in the excitation spectrum, in contrast to the I-phase, where the low-lying 
excitations are gapless. 

The systematic mean-field theory of the IC-transition was worked out by Bruce, 
Cowley and Murray in 1978 [2]. These authors distinguish IC-transition of type-I 
and type-II, depending on the absence or existence, respectively, of an inversion 
symmetry around g/p in the (disordered phase) soft mode dispersion. In the most 
simple case of type-I transition, condensation takes place only on wave vectors 
±qo For temperatures sufficiently below the mean-field transition temperature 
T c0 the system can be described, ignoring amplitude fluctuations, by the sinus- 
Gordon Hamiltonian. 

ft = 7 J d D x^(V<j ) -d)-g 2 vcosp4>} , (1) 

where 0(x) describes the long-wavelength distortions of the charge (spin, mass, 
flux etc.) density 

p(x)=p + Pl Re(e<^( x ))) . (2) 



Minimization of (0) yields for 5 < 5 C ~ gv l l 2 the solution of the C-phase (p = 
2n -n in = 0, 1, . . . ,p — 1). The solution for 5 > 5 C is a regular lattice of phase- 



solitons of distance I - , 
describe the I-phase close to the 
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and internal widths £o ~ ^/{.V9 V ) 1 > which 
C-transition. Far from the transition on has 
0(x) ~ <5x. Type-I transitions are therefore continuous. In general, a large 
number of C-phases is possible which may lead in certain lattice models to a 
devil's staircase behaviour of the modulation vector as a function of the misfit S 

[1]- 

Type-II transitions, on the contrary, are discontinuous and show in the I-phase 
an almost sinusoidal modulation of the order parameter. A recent example is the 
spin-Peierls system CuGeOs [3]. 

Thermal fluctuations have a strong effect on the IC-phase diagram in D = 2 
dimensions. If we exclude topological defects (i.e. vortices), then 

(i) the IC transition becomes an inverted Beresinskii-Kosterlitz-Thouless- 
transition with a reduced transition temperature T c (5 = 0) ~ Sirj/p 2 . 

(ii) Inside the I-phase solitons interact now by entropic repulsion. Close to the 
IC-transition the soliton-distance I increases as power law I ~ (5 — 5 C ) _ ^ S where 
Ps — 2(i-c) • H ere C is the thermal roughness exponent of a single soliton (domain 
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wall) Cth = (3 - D)/2 = 1/2. 

(iii) The spatial variation of the density 5p(x) = p(x) — p Q shows in the I-phase 
only quasi long range order (LRO): 

K(x) = ( e <^ x )-^(°))) ~ \x\- r >cos(27rz/pl) , (3) 

where S = 5e z . rj depends on I, T/7 and p and approaches 2/p 2 for I — > 00. 

(iv) Topological defects diminish further the ordered (C and I) phases and even 
change the topology of the phase diagramm for p < 4. In particular, for p < 2 
there is no direct IC-transition, both phases are separated by a fluid phase [4]. 
Closely results are expected for 1— dimensional quantum systems. It is interesting 
to remark, that a qualitatively similar picture emerges also in lower dimensions 
1 < D < 2, but with different singularities at the transitions [5]. 

On the contrary, in three (and higher) dimensions, thermal fluctuations have only 
minor effects on the phase diagram and are relevant essentially only in the critical 
region [6]. 

In the rest of the paper we investigate the influence of randomly distributed frozen 
impurities on the lock-in transition. Since impurities favour certain values of the 
phases, the impurity Hamiltonian can be written in the form 

H imp = 7/ d D xV(x,0) 
V(x,4>) = v / Acos(0(x) -a(x)), A = p\V*n im> (4) 

a(x) is a randomly frozen phase (0 < a < 2tt), 7V0 and n imp denote the strength 
and the concentration, respectively, of the impurities. The modell defined by ([J) 
and (gj) describes also an Xy-model in a crystalline (corresponding to a p-fold 
axis) and a random field. 

We will first discuss the case of a vanishing lattice potential, v = 0, and ex- 
clude topological defects for most parts of the rest of the paper. Larkin and 
in the present context first Fukuyama and Lee [7] have shown, that the im- 
purities destroy the translational LRO of the charge density wave on scales 
L ^> w A 1/(4_D) in all dimensions D < 4. Later studies[8] have shown, 
that for 2 < D < 4 K(x) decays as a power 



K(x) ~ e 



(5) 



with a universal exponent 77 = (^) e in lowest order in e = (4 — D). Thus, we 
regain now quasi-LRO in the I-phase in all dimensions 2 < D < 4. 
In systems in which the modulation of the I-phase has hexagonal symmetry, like 
in flux line lattices, the situation is more complicated. If we describe the distor- 
tions of flux lines by the displacement field u(x), the relevant correlation function 
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to describe long range translational order is given by ifc(x) = (e 4G ( u< - x * )_u< - - ) ) ), 
where G denotes the reciprocal lattice vector of the Abrikosov lattice. Recently, 
Emig at al. [9] found from a functional renormalization group (FRG) calculation 



K< 



LA*\ + Kzf) 5?lkj(x^ + zf) 



VG 

(6) 



where x = (x±,z), zi = {cu/ c^Y^z and k = c Ge /cu. c\\,c^ and cqq are the 
effective elastic constants of the flux line lattice (renormalized by thermal and 
disorder effects). In the marginal cases k = and k — 1 one finds from (|BD for 
the structure factor 

5(G + q) = / d 3 * e^ G (x) ~ ( q* + % 2 ) , (7) 

i.e. the structure factor exibits Bragg-peaks. The exponent 7/g is non-universal 
and depends on the value of k (1.143 < f)G < 1.159 for 1 > k > 0). In a large 
range of external fields k ~ 0o/16vrA 2 i? such that one could in principle test these 
predictions by measuring the field-dependence of width of the Bragg peaks. 
Next we come back to our original model (d), (P keeping the lattice potential 
v finite. Neglecting the non-Gaussian character of 0(x), which is justified for 
e <C 1, lowest order perturbation theory in v yields an effective Hamiltonian with 
a mass 



g 2 p 2 v <cosW>>~e- p2<< ^ 2> / 2 ~ f-^J . (8) 



-p 2 v/2 



Comparing this power of L with the L~ 2 behaviour of elastic term, we conclude, 
that the periodic perturbation is always relevant if p < p c = In this case, 
even an arbitrarily weak periodic potential will be relevant and we regain true 
translational LRO of the C-phase. 

For p > p c , on the other hand, weak periodic pinning is irrelevant, but for p 
close to p c we expect a transition to a commensurate phase for sufficiently strong 
v. A simple estimate for the threshold value v c follows from a comparison of 
the forces resulting from impurity and lattice pinning. With f imp « jL^ 2 and 
/„ m r yg 2 pv we get for the transition line v c = v c (A) ~ l/(g 2 pL 2 ) (or, inversely, 
A c = Ac(v)pa(gipv)^- D y 2 ). 

A more accurate description of the transition can be obtained from a functional 
renormalization group treatment [10], which confirms this estimate. The transi- 
tion turns out to be second order with a divergent correlation length £ ~ (v— v c )~ u 
in the C-phase, here v~ x = 4 — l)- Moreover, the order parameter for trans- 
lational LRO 

(1>) = {e-«)~(v-v e y > , P = ^ (9) 

is finite in this phase. On the contrary, in the I-phase the quasi-long range order 
of (B|) is regained. The fixed point, which describes this transition is at zero 
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temperature. The temperature eigenvalue — 6 = — 2 + e appears in the modified 
hyperscaling relation v{D — 6) = 2 — a typical for zero-temperature fixed points.. 




Figure 1: (a) Schematic RG-flow for 5 = and p > p c in the v — A plane. I,C 
and D denote the I-, C- and the disordered phase. (b)Possible phase diagram for 
finite 5 at fixed v for p > p c and (c) for p < p c 

As a site remark we note, that these exponents describe for p = 2 also the 
transition between the low temperature phase of the random field Ising model 
and the quasi-long range ordered phase of the random field XY -model. 
Apart from the change in the correlation functions at the transition there is also 
a change in the response on a small external drive / ex <C /imp, IV The creep 
velocity u creep can be written in the form 

Ucreepifc) ~ e~^7tT . (10) 

Here E c « 7^~ 2 , f c ~ max (/i mp , 7£~ 2 ) and fi = (D — 2)/2 and /i = — 1 for the 
I- and the C-phase, respectively. Because of the pronounced difference in the 
creep exponent /i in both phases measurement of the creep should give a clear 
indication about which phase is present. 

Recently measured I-V curves of the conductor o-TaS3 at temperatures below 
IK can be fitted by (10) with k = 1.5-2 [11]. The experimentally observed 
tendency to larger k for purer crystals confirms the above interpretation. In 
several materials, such as K0.3M0O3, the periods are near p = 4 commensurability 
at low temperatures. For this material one obtains £0 ~ 10~ 6 cm [10]. The typical 
parameters 7V0 ~ 10~ 2 eV, p\ = 10~ 2 and vp = 10 7 cm/sec yield, after proper 
rescaling of anisotropy, the estimate « 10 _4 cm for an impurity density of 
lOOppm. Thus it should in principle be possible to see commensurability effects 
if the misfit 5 becomes small enough, i.e. at low temperatures. 

So far, we have excluded topological defects. These can be considered if we treat 
0(x) as a multivalued field which may jump by multiples of 2n along certain 
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surfaces. These surfaces are bounded by vortex lines. For 5 — = v it has been 
shown recently that for weak enough disorder strength A < Ap, the system is 
stable with respect to the formation of vortices [12]. However, vortex lines will 
proliferate for A > A^. At present it is not clear whether the corresponding 
transition is continuous or first order. For 5 = but v > we expect that 
this transition extends to a line A D (v) until v reaches a critical value vd with 
Az)(t>£)) = A c (vd) (see Figure 1). For larger v the transition is probably in the 
universality class of the p-state clock model in a random field, which has an upper 
critical dimension D c = 6 . A non-zero value of 5 will in general increase the size 
of the incommensurate phase, as schematically scetched in Figure 1. 
Our results can also be applied to the pinning of flux line lattices in type-II 
superconductors. In layered superconductors, the Cu02 planes provide a strong 
pinning potential favoring, for flux lines oriented parallel to the layers, a smectic 
phase with translational order present only along the layering axis [13]. The 
influence of disorder on this phase is described by the CDW model studied in 
this paper, if the CDW phase 0(x) is identified with the deviations of the smectic 
layers from their locked-in state. Also a FLL oriented perpendicular to the layers 
in general feels a weak periodic potential of the underlying crystal, but now the 
flux line displacements are described by a vector field. Since also in this case weak 
disorder leads only to logarithmically growing transverse displacement of the flux 
lines [9], a disorder driven roughening transition of the CDW type studied above 
can be expected for the flux line lattice. 
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